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A conjectured connection to quantum gravity has led to a renewed interest in highly damped 
black hole quasinormal modes (QNMs). In this paper we present simple derivations (based on the 
WKB approximation) of conditions that determine the asymptotic QNMs for both Schwarzschild 
and Reissner-Nordstrom black holes. This confirms recent results obtained by Motl and Neitzke, 
but our analysis fills several gaps left by their discussion. We study the Reissner-Nordstrom results 
in some detail, and show that, in contrast to the asymptotic QNMs of a Schwarzschild black hole, 
the Reissner-Nordstrom QNMs are typically not periodic in the imaginary part of the frequency. 
This leads to the charged black hole having peculiar properties which complicate an interpretation 
of the results. 
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A. The quasinormal modes 

Black holes oscillate. The associated quasinormal modes (QNMs) of oscillation are relevant for many reasons. Most 
^ importantly, numerical relativity has provided ample evidence that the QNMs dominate the gravitational- wave signal 
associated with many processes involving dynamical black holes (such as the formation of black holes in gravitational 
7-H , collapse or binary merger) . Since the QNMs encode information concerning the parameters of the black hole one may 
J> ' hope that the gravitational-wave detectors that are now coming into operation will be able to use these signals to 
investigate the black hole population of the Universe. 

Even though most studies of QNMs have been motivated by their potential astrophysical relevance, there are 
several other reasons why one might be interested in understanding the spectrum of oscillations of a black hole. 
In particular, the modes have played a key role in discussions of black hole stability Jj. A closely related issue 
concerns mode-completeness. It now seems clear that the QNMs do not form a complete set (at least not in the 
I conventional sense) because of the presence of power-law tails caused by the backscattering of waves (see |2| for 
a discussion and references). The problem also has interesting computational aspects. While the slowly damped 
QNMs, for which |Re bjM\ » |Im ijjM\, are relatively easy to compute, highly damped modes present a challenge. 
I • The main difficulty concerns the fact that, in the frequency domain, the QNM eigenfunctions (which represent purely 
outgoing waves at spatial infinity and purely ingoing waves crossing the event horizon) grow exponentially. This 
• • means that one must, in principle, achieve exponential precision in order to impose the boundary conditions. The 
' first reliable calculation of high QNM overtones was performed by Leaver in the mid-1980s using continued fractions 
, [^]. An alternative approach to the problem proceeds via analytic continuation using complex coordinates, borrowing 
^ ' standard asymptotic techniques from quantum mechanics. In particular, this was the fundamental idea behind the 
5^ , numerical phase-amplitude method that has been used to calculate high precision QNMs of both Schwarzschild and 
Reissner-Nordstrom black holes |^^. Both these results, and those of NoUert and Schmidt |^, agreed with those of 
Leaver. 

While the reliability of Leaver's method is now well established, it was not without controversy a decade ago. The 
discussion concerned the behaviour of the QNMs in the limit of very high damping. Leaver's results (for the first fifty 
modes of a Schwarzschild black hole) indicated that the modes would asymptotically behave as 

w„M ~ 0.08- ^(2n-h 1) asn^oo 

This was contradicted by results obtained by Guinn et al using a WKB formula ^ . Their calculation suggested that 
the real part of the QNM frequencies would vanish asymptotically. The controversy was resolved by two calculations 
which agreed that the correct asymptotic result was 

- 0.04371235- ^(2n+ 1) as n ^ 00 (1) 

The first calculation was based on a slight reformulation of the continued fraction algorithm Q , while the second used 
a high order phase-integral formula j^, il^] . The latter calculation also shed light on the reasons for the breakdown of 
the method used by Guinn et al. 0- This issue is further discussed in Refs. pdllT^ . 
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B. Is there a quantum connection? 



Despite the fact that the laws of black-hole thermodynamics are by now well established, many issues remain 
unclear. It is, for example, not clear to what extent the black hole entropy 

^-I (2) 

where A is the area of the event horizon, can be understood in terms of the statistics of a given set of microstates. 
(We use units such that c — G = h = 1 throughout the paper.) Bekenstein and colleagues have discussed this problem 
in terms of a quantised area (see for example llSl') . Then, in analogy with a typical finite system, a black hole would 
have a discrete spectrum. One can argue that 

A^nAlnk with 7i 1, 2, 3, ... and fc = 2, 3, 4, ... (3) 

Comparing this with 

A = 16nM^ -> AA = 32ttMAM = 327rMw (4) 

where we have associated the "energy spacing" with a frequency through AM = AE = uj (roughly speaking, uj 
corresponds to a "transition energy" ) , one finds that the spacing between consecutive states (for macroscopic black 
holes, with M >> fi) will correspond to a frequency 

UJ — (5) 

The standard argument, which favours k ~ 2, would make the entropy spacing between energy levels exactly one 
"bit", which is attractive from the information theory point of view. A key result is that, in this picture any radiation 
will be emitted in multiples of the fundamental frequency u. Hence, essentially no radiation should be radiated with 
frequencies below w. If true, this is a highly significant conclusion since it suggests that the quantum nature of black 
holes might be observable at the macroscopic level. 

As has been demonstrated by Ashtekar and his collaborators one can arrive at the same conclusions within 
the framework of loop quantum gravity (see [T3 | for a nice introduction). This approach is based on the notion of 
quantum geometry, which means that it is natural to ask what the quantum of area might be. Again, it is possible 
to draw conclusions from considerations of macroscopic black holes. For large black holes it has been shown that 0| 

A = n7B/47r\/3 (6) 

The parameter jbj is an unknown "natural constant" called the Barbero-Immirzi parameter. It plays an important role 
because it fixes an ambiguity in the theory. If this parameter could be determined by an independent "experiment" the 
theory would become predictive. In fact, the Barbero-Immirzi parameter can be fixed by carrying out a calculation of 
the black-hole entropy and comparing to the standard result, Eq. (0). By attributing entropy to microstates (actually 
nodes of spin networks, see |15j) it can be shown that 

ln2 ^ ln2 , . 

-A 7s/ = — y= (7) 



provided that the lowest permissible spin is 1/2. Again the area quantum is AA = 41n2. 

So what does this have to do with the black hole QNMs? The possible connection between the classical vibrations of 
a black-hole spacetime and various quantum aspects has been discussed for quite some time. For an early contribution 
to the debate, see Ref. 0|. The recent interest in a possible association between the two problems followed a very 
simple observation. A few years ago. Hod 0| noticed that the numerical results for asymptotic Schwarzschild QNMs 
seemed to suggest that 

In 3 

Re LuM as |Im cj| ^ oo (8) 

If we identify the (real part of the) asymptotic QNM frequency with the quantum interstate spacing, we can use this 
value in to get 

u; = ^ A^ = 41n3 (9) 

8nM ^ ' 
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This is a tantahzing result. It would fit nicely into Bekenstein's thermodynamical pict ure provided that the "fun- 
damental" frequency is associated with fc = 3 rather than 2. Furthermore, Dreyer [l^ has shown that this would 
be the prediction of loop quantum gravity if it were based on SO (3) rather than SU(2). In this case the predicted 
Barbero-Immirzi parameter would be 

In 3 , , 

There are, of course, problems associated with this change. In particular, SO (3) is not favoured because it would seem 
not to allow coupling to fcrmions. Hence, one would need to either explain why fermions should be excluded |l9j |. or 
provide an alternative derivation of the black hole entropy (perhaps using a different statistics for the quantised area 
states ^3)- Various other relevant issues have been discussed in Refs. |2lll23 .l2 ^ . 

On the other hand, we now have a testable prediction. If Hod's argument holds, one should be able to learn 
something useful from, for example, the asymptotic behaviour of the Reissner-Nordstrom QNM frequencies. However, 
until very recently only the slowly damped QNMs of Reissner-Nordstrom black holes had been calculated |^ 12^. I25L 

El 113. 

The present paper is motivated by the recent discussion [2^ |2^ [s^, |^ H^l- Our aim is to use the complex 
coordinate WKB method |U |3 |33 in its very simplest form to determine asymptotic QNMs for Schwarzschild 
and Reissner-Nordstrom black holes. As we will see, this leads to results that agree with those of Motl and Neitzke 
[2^. Furthermore, our analysis fills several gaps left by their study and extends the discussion of the final result 
considerably. 



11. KEY PRINCIPLES OF THE COMPLEX COORDINATE WKB ANALYSIS 



The equations governing various classes of non-rotating black-hole perturbations can be written |3 

dr'i 



f^ + [u;'-V{r)]i;^0 (11) 



where we have assumed that the perturbations depend of time as e The tortoise coordinate is defined by 

where A depends on the spacetime geometry, and is here given by 

A = - 2Mr + (13) 

where M is the mass of the black hole and q is its electric charge. The two solutions to A = determine the location of 
the event horizon = M + \J — and, for charged black holes, the inner Cauchy horizon r_ — M — Af ^ — q^ . 
The definition H12|) is such that the causally attainable spacetime region outside the black hole, r_|_ < r < cx), is 
mapped onto —00 < r» < 00. The effective potential V{r) is of short range, which means that V ~ 6=^*"'^* both near 
the horizon and at infinity. With our chosen time-dependence the solution behaving as e*"''* represents an outgoing 
wave at infinity, while e~^^'^* is a ingoing wave near the horizon. 

It is easy to explain one of the main difficulties associated with the QNM problem. Suppose we want to calculate a 
damped QNM, i.e. a solution for which Im loM < 0. Then the solution that represents outgoing waves at infinity will 
grow exponentially as r, increases. This means that we will need exponential precision in order to filter out the ingoing- 
wave contribution and impose the desired boundary condition. This is not a straightforward computational task. A 
similar difficulty arises with the boundary condition at the horizon. However, the problem becomes straightforward 
if we analytically continue into the complex coordinate plane. Suppose, for example, that we analyse the problem 
along a line such that tur^ is purely real. Then the two asymptotic solutions would be purely oscillatory and it would 
be easy to impose the boundary conditions. Of course, before we can benefit from this idea we must understand 
how the solutions change under analytic continuation. Fortunately, the relevant principles are well known from 
WKB /phase- integral theory. In fact, in this paper we will only use results that were well understood at least 40 years 
ago 32., ,3^4,, ,3,5] . The advantage of this kind of analysis is that it does not require the use of complicated comparison 
equations to spot special solutions whose analytic properties can be exploited (it is an "atomic" description rather 
than a "molecular" one). 
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In order to analyse the black hole problem we prefer to work in the complex r-plane (this is natural since the 
tortoise coordinate is multi- valued) . Introducing a new dependent variable 



we readily rewrite the perturbation equation as 



A 



1/2 



i?(r)* ^ 



where 



d_ / A 

dr \ 



1 A^ /A~ 

2 dr'^ \ 



As is well known, the two WKB solutions to an equation of form (|15|l can be written 



/SW = Q~'/'Mexp 



±i 



Q{r')dr 



(14) 



(15) 



(16) 



(17) 



with = R. However, under some circumstances it is useful to use a slightly different function as basis for the 
approximation. As we will discuss in the next section, we will exercise this freedom in our analysis of the asymptotic 
black hole problem. Note that, without loss of generality, the lower limit of integration t is customarily taken to be 
one of the zeros of Q. Throughout this paper we will indicate the relevant lower limit of integration by a subscript 
on /i_2 as in ((T7jl . 

The zeros and poles of the function Q play a central role in any complex coordinate analysis of l|15f) . From each 
simple zero of emanates three so-called "Stokes lines". Along each of these contours Q{r)dr is purely imaginary, 
which means that one of the two solutions grows exponentially, while the second solution decays, as we move away 
from t. In other words, one of the solutions is exponentially dominant on the Stokes line, while the other solution 
is sub-dominant. Analogously, one can define three "anti-Stokes lines" associated with each simple zero of . On 
anti-Stokes lines Q{r)dr is purely real, which means that the two solutions are purely oscillatory. As we cross an 
anti-Stokes line, the dominancy of the two functions /i,2 changes. The three-fold symmetry associated with each zero 
of is clear from Figs. Q and El 

Stokes lines are vital for WKB analysis, because it is in the vicinity of these contours that the solution changes 
character. That is, if the solution is appropriately represented by a certain linear combination of /i and /2 in some 
region of the complex r-plane, the linear combination will change as the solution is extended across a Stokes line. The 
induced change is not complicated: The coefficient of the dominant solution remains unchanged, while the coefficient 
of the solution which is subdominant on the relevant Stokes line picks up a contribution proportional to the coefficient 
of the dominant solution. This is known as the "Stokes phenomenon" [37'|. The contant of proportionality is known as 
a "Stokes constant" . This change is necessary for the particular representation H17|) to preserve the monodromy of the 
global solution. Terms that are exponentially small in one sector of the complex plane may be overlooked. However, 
in other sectors they can grow to exponentially dominate the solution. By incorporating the Stokes phenomenon, we 
have a formally exact procedure which leads to a proper account of all exponentially small terms. 

In the particular case of an isolated simple zero t of the problem is straightforward [ssL l34l |. Suppose that the 
solution in the initial region of the complex plane is given by 



Then, after crossing a Stokes line emanating from t (and on which /i is dominant) the solution becomes 



* = c/f ^ ± icf, 



(t) 

2 



(18) 



(19) 



The sign depends on whether one crosses the Stokes line in the positive (anti-clockwise) or negative (clockwise) 
direction. It is crucial to note that this simple result, i.e. that the Stokes constant is ±i, only holds when the Stokes 
line emanates from the zero that is used as lower limit for the phase-integral. That is, when we want to use the above 
result to construct an approximate solution valid in various regions of the complex plane we must often change the 
reference point for the phase-integral. This leads to the need to evaluate integrals of the type 



1^3 



Q{r)dr 



(20) 
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where ti and tj are two simple zeros of Q^. 

A final issue that must be mentioned before we proceed concerns branch cuts. In general we need to introduce 
branch cuts from the simple zeros of in order to ensure that the phase integrands remain single-valued. However, 
one can usually place these cuts in such a way that they do not affect the analysis. In the following derivations this 
is the case. We will choose the phase of the square-root of such that 

Q ^ i?i/2 _ ^ asr^oo (21) 

This means that the outgoing-wave solution at infinity is proportional to /i while the ingoing-wave solution at the 
horizon is proportional to f2- 



III. THE SCHWARZSCHILD PROBLEM 



A. Evaluating the phase-integrals 



In the case of Schwarzschild black holes (when q ~ 0) there are two classes of gravitational perturbations, usually 
refered to as axial and polar [3^ . We will only consider the axial case here. This is, however, no restriction since it 
has been shown that the two cases are isospectral 36j. That is, the QNMs are the same in both cases. Once we have 
derived the relevant WKB condition for axial gravitational perturbations we will discuss the case of electromagnetic 
waves in the Schwarzschild background. 

In the case of axial gravitational perturbations, the effective potential is 



{1(1 + 1) 



6M 



Z = 2,3,4, 



(22) 



From this and Eq. we deduce that the function R has two second order poles and four zeros. Since the zeros are 
closely associated with the Stokes phenomenon, we need to know their location, as well as the nature of the Stokes 
and anti-Stokes lines. It is easy to show that, when Im co — > —oo the zeros all approach the origin of the complex 
r-plane. This allows us to simplify the analysis considerably. Expanding in a power series near r = we have 



R. 



4M2 



IBM 



21 



15 

4^2 



(23) 



sufficiently near to the origin. Note that this approximation contains no reference to I. As — s- oo the l-dependent 
terms are only higher order corrections. 

Given this behaviour it is easy to show that the exact solutions to IjlSfl should behave like 



r^/'^'^^ as r 







(24) 



Meanwhile, if we take — R we get 



-1/2 

Qdr 



^-1/2 

±t\ — mr 
V 4 



(25) 
(26) 



which means that 



/i 



(27) 



In other words, the approximate solutions do not have the correct behaviour in the vicinity of the origin. However, 
by choosing 



1 

¥2 



(28) 



we obtain approximate solutions with the desired behaviour near the origin. This is analogous to the "Langer 
modification" l{l + 1) ^ {I + 1/2)^ that is used in the WKB analysis of radial quantum problems. In principle, one 
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can also adjust the approximation near the poles at r± in the black hole problem, see for example but since we 
are assuming that \lu\ is large such alterations would not affect our final result. Hence, we will use 



4M2 



to 



(29) 



In Fig. n we show the anti-Stokes and Stokes line geometry pertaining to (|28() for large The zeros have been 
labelled in the same way as in Refs. 0,^2- o^'' analysis of the QNM problem we will need the phase-integrals 713 
and 723. That is, we need to evaluate 



/ = QdrK± 



r 
2M 



21 1/2 



dr 



(30) 



where the limits are two neighbouring zeros of Q^. Letting y — ojr^ /AM the zeros map to —1 or 1 and this integral 
becomes 



/ = ± 





1 ■ 








dr^± I 




y . 


'J in 



— — - — dx = ± [sinhx — 2 arctane^]°_ = ^tt 
cosh X 



(31) 



Hence we see that, up to the sign, the integrals we need are identical. Furthermore, it is easy to show that with our 
chosen phase for Q we obtain 



7 ^ 



-713 



-723 = 



(32) 



B. A WKB condition for asymptotic Schwarzschild QNMs 



We now combine the monodromy argument of Motl and Neitzke j23| with the standard complex coordinate WKB 
results described in the previous section. This will provide a clear argument in support of the known result for 
asymptotic Schwarzchild black hole QNM frequencies. 

For frequencies such that |Im uj\ >> |Re iu\ the pattern of Stokes and anti-Stokes lines for the Schwarzschild problem 
is as sketched in Figure^ Assuming that Re loM > the outgoing wave boundary condition at spatial infinity can be 
analytically continued to the anti-Stokes line labelled a in the figure. This issue is discussed in detail in, for example, 
[Tll |. In order to obtain a quantisation condition for highly damped QNMs, we analytically continue the solution along 
a closed path encircling the pole at the event horizon. This contour starts out at a, proceeds along anti-Stokes lines 
and account for the Stokes phenomenon associated with the zeros ti, t2 and ^3, and eventually ends up at a. In the 
analysis we will assume that all zeros and poles of are isolated and can be accounted for individually. 

With the chosen phase of Q the outgoing-wave solution at point a is 

V'a = /^^ (33) 

Since no Stokes lines cross this contour, this solution will not change in character along the anti-Stokes line that 
connects point a with the zero ti. This means that we can readily extend the solution to the vicinity of ti. However, 
if we want to extend the solution to point b on a neighbouring anti-Stokes line we must account for the Stokes 
phenomenon. With our choice of phase for Q, the function f^^^ is dominant on the Stokes line which we must cross 
in going from a to h. This means that we will get 

i^b = /^^ - (34) 

since the Stokes constant for a single well separated zero is —i if we move around the zero in the clockwise direction. 
Changing the lower limit of integration to ^3 we get 

^f, = e''^"/f ^'^ - ie-''^"/i*'^ = -.ff + if 2^^ (35) 

Now extending this solution to point c we do not cross any anti-Stokes lines so /i remains dominant. Hence, at c we 
obtain 

V'c = -.f^^ + 2*/^) (36) 
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FIG. 1: A schematic illustration of the Stokes (dashed) and anti-Stokes (solid) lines for the Schwarzschild problem in the 
complex r-plane. We have assumed that the frequency is large and purely imaginary. For frequencies with a small real part, 
the pattern of Stokes and anti-Stokes lines changes only slightly (cf. various figures in Hill). The open circles represent the 
four zeros of , while the filled circles are the two poles (at the origin and the event horizon, respectively). The outgoing- wave 
boundary condition at infinity is imposed on the anti-Stokes line labelled a. For uj = a — if3, with [a, /3] > this anti-Stokes 
line asymptotes to a straight line at an angle arctan/J/a, i.e. in the first quadrant of the complex r-plane. 



Since we do not cross any Stokes lines in moving from c to c', the linear combination remains a valid solution. 
However, we need to note that the phase-integral is now evaluated along a contour that loops around the pole at 
r = 2M. If we replace this integration contour by one that lies to the left of the pole, we get 

V-c = -e^^/f^^ + 2ze-*7i*^^ (37) 

where F is the integral of Q along a contour encircling the pole at r = 2M clockwise. 

We now want to connect the solution to the point d. In order to do this we must first ensure that the lower limit 
of the phase-integrals is ^2- That is, we use 

V-c' = -e*(^+'^^-^Vf + 2te-'^^-'^'^^ f^''^ = e'^/f - 2te-'^ f^''^ (38) 

When we step inside the large anti-Stokes lobe in Fig. ^ the dominance of /i and /2 is interchanged. Then, crossing 
the first of the two Stokes lines inside this lobe we obtain 

= - [e^^ + 2e-^n fi"^ + 2^e-'^ft^ (39) 

Connecting this solution back to c, again crossing a Stokes line where /i is subdominant, we get (using a bar on c 
to denote the fact that we have encircled the pole at the origin) 

= - [e^r + 4g-.r] ^(^3) ^ 2^e-^^(*^) (40) 
Now reversing our steps and connecting this solution to b we get 

^ _ [e«r ^ 4g-.r] fits) _ ^ ^^r ^ ^^-.r^ e'^^ f^^ 

= [e'^ + 4e-'^n /^^ + i [e^i^ + 26"^^ jiti) (4^) 
and (finally) returning to a we have 

= [e^r ^ 4g-.r] fit.) ^ 2, [e^r ^ g^-.r] fit.) (42) 
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Comparing this solution to H33|l we clearly must require that 

+ ie~'^ = (43) 
in order for the two solutions to be the same. This would lead to 

^-a = e-^r^^) (44) 

from which we see that the clockwise monodromy of this solution is e~^^ . 

Let us now perform the same analysis in the vicinity of the pole bX r — 2M. With our choice of phase for Q, the 
solution that represents "ingoing waves" near the event horizon is 

i^H = Pft^ (45) 

For this solution it is trivial to show that the clockwise monodromy is (again) e^*^. 

Since a necessary condition for the two solutions H33(l and H45|l to represent the same global solution to the problem 
is that they have the same monodromy, we conclude that H43(l is the appropriate WKB condition for highly damped 
QNM solutions. 

Finally, using the fact that 

F = (l Qdr = -27rz Res Q = -A-rriuM (46) 

/ r=2M 

our WKB condition can be written 

e'^'^ = e^""^^'' = -3 (47) 

and it immediately follows that 

uoM = — ln3 — -(riH — ) asn^oo (48) 



?7r 4 V 2, 

This is the desired final answer, in complete agreement with l30l |. Of course, it is worth noting that our 

derivation is conceptually very simple as it only appeals to basic WKB principles. 

It is interesting to discuss other classes of black hole perturbations. In order to do this we note that, had we not 
used the particular value 7 = tt in our derivation, we would have arrived at the condition 

gSTT^M ^ _ 2cos27 (49) 

We can easily use this condition to determine the asymptotic QNMs for electromagnetic waves. 

In the case of electromagnetic waves propagating in the Schwarzschild geometry the relevant effective potential is 

V = l{l + \)^ (50) 

From this we find, using iQlI) and idHl, that 

4M2" 



4M2 

The topology of the problem is still represented by Figure and given H31() we find that 

1/2 

t,r 

2M 



(51) 



Qdr^± I ^ 



dr = ±1 (52) 



Hence, we have 7 = 7r/2 and it follows immediately from H49|) that the real part of the QNM frequencies vanishes 
asymptotically. 
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IV. THE REISSNER-NORDSTROM PROBLEM 



Having verified the known result for Schwarzschild black holes we will now consider the QNMs of the Reissner- 
Nordstrom geometry. This problem has recently been discussed by several authors, see [I^I^UjI, but the implications 
of the results still seem far from clear. The Reissner-Nordstrom problem is more complicated than the Schwarzschild 
one because of the presence of the inner horizon as well as two additional zeros of the relevant . Nevertheless, the 
analysis proceeds almost exactly as in the previous section. 

In the case of charged black holes, one still only needs to consider axial perturbations. Just as in the Schwarzschild 
case the axial and polar perturbations are isospectral For axial perturbations one has two distinct effective 

potentials, cf. p^. 



(53) 



where 



/3i^2 = 3Af T W + 4(Z -!)(; + 2)(?2]i/2 



(54) 



In the Schwarzschild limit (g 0) these potentials approach pure electromagnetic (Vi) and gravitational (V2) pertur- 
bations, respectively. In the general charged case, the solution corresponds to coupled electromagnetic and gravita- 
tional waves. That this should be the case is natural since oscillations of a charged gravitational field will inevitably 
generate electromagnetic waves, and vice versa. Note that the g ^ limit is, in fact, singular due to the coalescences 
of poles and zeros that change the character of the function Q^, cf. Figs. ^E^ndl^l 



A. Evaluating the phase-integrals 



When <7 7^ the function R has three second order poles and six zeros regardless of whether we consider V\ or 
V2. Furthermore, just as in the Schwarzschild case, the zeros all approach the origin of the complex r-plane when 
\u>\ — > 00. Again expanding in a power series we find 



6?" 



6 



(55) 



in the immediate neighbourhood of r = 0. The leading order behaviour is the same for both classes of perturbations 
(Vi and V2). Hence, our analysis will hold for both sets of perturbations and the final result will be identical in the 
two cases. 

Repeating the argument from the previous section one can show that the choice (|28|) still leads to the WKB solutions 
having the behaviour expected of the exact solutions near r = 0. Thus wc will use 



25q* 
4^6 



(56) 



As in the previous section we will need the phase-integral connecting neighbouring zeros of Q^. That is, we require 

1/2 



1= Qdr^± - 



25g^ 



4^6 



dr 



(57) 



Letting y — 2ujr^/5q^, and using two neighbouring zeros of as limits, the integral becomes identical to that of the 
Schwarzschild problem (apart from a multiplicative factor), and we find 



J = ±- 



57r 



(58) 



Yet again all the integrals we need are the same (up to sign). Chosing the phase of Q as in the previous case, and 
labelling the zeros as in Fig. |21 we have 



-712 



~732 = 743 



-754 



5tt 
~6 



(59) 
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Let us denote the integral along a contour that encircles the pole at r+ (in the negative direction) by Fg. Then 

Te = 725 + 754 + 743 + 732 = 725 + 743 - "^1 (60) 

where a tilde indicates that the integral is taken along an anti-Stokes lobe to the right of either the pole at the event 
horizon (725) or the pole at the inner Cauchy horizon (743), cf. Fig|21 Similarly, we define 

Ti = -743 + 743 = -743 + 7 (61) 

These definitions allow us to write 

743 = -r, + 7 (62) 

and 

725 = Te + r, + 7 (63) 
The two integrals Fg and F^ are readily evaluated using the residue theorem. We find 

F- = -O.TTi Rp.R Q = -7ri + = -iri (l + nf (64) 

where we have defined the dimcnsionless parameter 



g2 



M' 

The integral around the inner horizon is 



'^=\/l-172 (65) 



Vi = —2TTi Res Q = iri — - = ni (1 — k)^ (66) 



r— r_ 



B. The WKB condition for Reissner-Nordstrom QNMs 




FIG. 2: Stokes and anti-Stokes lines for Reissner-Nordstrom problem. Just as Fig. this is a schematic picture based on a 
purely imaginary frequency. The three poles in the problem are represented by filled circles and correspond to, from left to 
right, the origin, the inner horizon and the event horizon. The six zeros of are shown as open circles. 
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The pattern of Stokes and anti-Stokes lines for the Reissner-Nordstrom problem for |Im u>\ » |Re is shown in 
Figure|21 Introducing the outgoing- wave boundary condition on the appropriate anti-Stokes line (at point a in Fig.|2J), 
just as in the Schwarzschild problem, and then moving around the pole at r = r+, taking full account of all involved 
Stokes phenomena (assuming the zeros are well isolated and can be treated individually), we arrive at the following 
WKB condition for highly damped QNMs 

e^^r^ = 1 - (1 (l + e^'^) (l + e-^^.) (57) 

(the complete derivation is provided in the Appendix). Since we know that 7 = Stt/G this becomes 

e^^^e = 1 _ 2 (^1 + cos y ^ (1 + e-^^^') = -2 - Se-^^r. (gg) 

which can be shown to be identical to the condition recently derived by Motl and Neitzke '2^ using matched asymp- 
totics. That this condition agrees well with the numerical solution to the QNM problem, for various given overtones, 
has been shown by Berti and Kokkotas j8ll |. 



C. Approaching the Schwarzschild limit 

It is interesting to consider what happens when the black hole approaches the Schwarzschild limit, since H68|l is 
clearly at variance with the result (|47|l for Schwarzschild black holes. For — s- cx) then g — s- 0, we have 

T,^-Amu:M(l + -^\ (69) 



and 



Thus (|68|l predicts that we ought to have 



16M4 



gSTT^A/ = -5 (71) 



Yet the Schwarzschild calculation predicted that the right-hand side should be —3! The reason for the discrepancy 
is, however, easy to explain. As we have already mentioned the Schwarzschild limit is singular. At our level of 
approximation one cannot move from the topology of Fig.|21to that of Fig. ^ in a non-trivial way. This would require 
a uniform approximation involving the coalescence of two zeros and two poles. 

The asymptotic behaviour of the QNMs of a charged black hole is always given by (|68|l . and hence corresponds to 
the real part of the QNM frequencies approaching IuS/Stt in the limit of infinite damping. But there is also likely 
to be an intermediate regime in which the highly damped QNMs more resemble the Schwarzschild result, i.e. where 
Re bjM w ln3/87r. As the Schwarzschild limit is approached, this latter regime tends to dominate, with the true 
Reissner-Nordstrom asymptotic behaviour being relevant only for extremely rapidly damped QNMs. That this is the 
case can be understood by the following argument: Let us consider a black hole with an infinitesimal charge and 
QNMs such that |Im loM\ >> |Re ljjM\. The main difference between this problem and the Schwarzschild case is the 
presence of the double pole associated with the inner horizon r_, and two additional zeroes of the function R, defined 
by (|16|) . As the imaginary part of the QNMs increases, all six zeros of R move towards the double pole at the origin. 
Our analysis of the problem is only relevant when the topology is that illustrated in Fig. |21 i.e. when the pole at r_ 
lies outside the circle on which the six zeroes are located. For infinitesimally charged black holes there will also exist a 
regime where the topology of the the four zeroes already present in the Schwarzschild case are essentially unchanged. 
For this to be true, we need to have 

1^1 « {M/uf^ (72) 

At the same time the pole at r_ , and the two additional zeroes that come into existence when the black hole attains 
charge (and which emerge from the origin together with r„), must lie well inside the circle of "Schwarzschild" zeroes. 
This corresponds to 

|r_| « \t\ ^ q2 « 2M{M/Lof'^ (73) 
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If these conditions are met one would expect the QNMs to be similar to the Schwarzschild ones. Figure|31is a schematic 
illustration of the two regimes for highly damped Reissner-Nordstrom QNMs. Even though it is easy to explain the 
behaviour in principle, it is not straightforward to extend our analysis into a consistent scheme for calculating the 
QNMs in the intermediate regime. The main reason for this is the need to evaluate the phase-integral 7. In the 
intermediate regime we can no longer make use of the power series expansion around the origin that led to H29() . Such 
an expansion is only valid up to the nearest pole, i.e. r_. Resolving this issue, and determining the QNMs also in 
this intermediate regime, may be an interesting problem but we do not consider it further here. 



D. Solving the QNM condition 

In order to discuss the solutions to the Reissner-Nordstrom QNM condition, it is useful to rewrite (|()8|l as 

^SttujA! 2g^27rci;M(l — k)^/k (74) 

where we have used ()64ll and 1)66(1 . From this we easily see that in the extremely charged case, as k — > 0, we regain 
the Schwarzschild result [s^- That is, we get 

Re LoM as g ^ M (75) 

It is, however, not clear to what extent this result is relevant. After all, the extreme Reissner-Nordstrom limit is 
singular in the sense that the two poles at r± coalesce as g — > M. This means that the topology is no longer that 
illustrated in Fig. El and hence the extreme Reissner-Nordstrom case would require a separate analysis. 
In order to analyse the general case, we introduce 

(1 - k)'^ 

y = SttcoM and k = ^ , ' (76) 

4:K 

The condition can then be written as 

e« = -3 - 2e-''y . (77) 

In general, this condition must be solved numerically. One way to do this is to first separate the real and imaginary 
parts of the equation. Letting y — a — i(3 we obtain the two equations 

e"cos/3 = -3-2e-'="cosfc/3 (78) 
e"sin/3 = 2e"''"sinfc/3 (79) 

These equations are very useful. First of all we see that, if we want the solution to be periodic in /3, we must require 
simultaneously 

P ^ l3 + 2mT and fc(/3 + 2mr) = kf3 + 2rmr (80) 
where n and m are integers. That is, periodicity in Im lo is only possible if 

k = (81) 
n 

Since < fc < cxd we see that we will pass through all (positive) rational values for k as the charge of the black hole 
is varied. Hence, there will be an infinite set of cases where the solution is periodic in the imaginary part. However, 
whenever k is not a rational number, the spectrum ceases to be periodic in the imaginary part. This is a significant 
observation because it illustrates that the asymptotic Reissner-Nordstrom spectrum is generally very different from 
that of a Schwarzschild black hole. One might wonder if the inclusion of higher order terms in the approximation 
would reimpose periodicity in Im to for irrational values of fc. However, this seems unlikely since it would require a 
surprising "conspiracy" between these terms. The difficulty (due to truncation error) of representing exactly rational 
numbers would pose significant challenges for a numerical verification of this result. 
For m and n integers we can introduce z — e^/" to get 

+ 3z" + 2 = (82) 
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which obviously has n + m roots (one of these polynomial cases was discussed in ^525). These roots lead to 

SttujM = y = nhi\z\ + in argz + 2np'Ki p = 0,l,2,... (83) 

Only some of these roots will be compatible with the derivation of our WKB condition since we assumed that 
Re ujAI > at the outset. The permissible roots correspond to the smallest basis set required to recontruct the 
spectrum. Repetitions of these roots on other sheets of z (different k), for large p generate the asymptotic spectrum, 
cf. figure|3| Hence, it is interesting to consider some of these polynomial solutions. In Table^we give the roots that 
are compatible with the derivation of our QNM condition for some selected values of fc = m/n. 



In 3 / 8ff In 5/ 



-0.625 




-D.625 



RcaiM 



Basis set of 
■ eigenvalues foL 
fixed Lational jt. 



Repetition of 
basis set of 
eigenvalues foL 
fixed lational jt 
at hishei lm6JM. 



-Im aJM 



SCdW 



FIG. 3: Left panel: A schematic illustration of the behaviour of the first few solutions to the Reissner-Nordstrom QNM 
condition Hfi8|l . As the charge is increased the QNM frequencies spiral towards InS/Svr — (2p + l)i/8. The figure illustrates 
how the spirals tighten with increasing damping. Generally, the asymptotic QNMs correspond to the high damping limit of 
the figure. Middle panel: In the case of fixed rational k = m/n, the spectrum is divided into repetitions of basis sets of 
eigenvalues from low-lying states (indicated by dots). The number of modes in the basis set increases with n. For irrational k 
this periodicity is destroyed. Right panel: A schematic illustration of the two asymptotic regimes of the Reissner-Nordstrom 
problem, RN and SCHW in the figure. In the limit of infinite damping the QNMs are determined by H68|l . but there also exists 
an intermediate "Schwarzschild-like" regime . The boundary between the two regimes is (roughly) given by « 2M (Af/ 



The data given in Table provide the basis needed to generate the asymptotic QNM spectrum in the simplest 
periodic cases. As shown in the table, there are n distinct QNMs in each basis set. These yield the asymptotic 
spectrum for large values of p, cf. H83|l . These results are only valid in the regime where the topology of the problem 
is that shown in Figure |21 From our discussion of the Schwarzschild limit above, we deduce that this is the case when 



|Im ujM\ » ^M/qf 



(84) 



It is relevant to ask if there are any values of k for which we have purely imaginary frequencies. One can show that 
a = is only compatible with (|78|) and H79|) provided that 



That is, we must have 



sin/3 = and cos/3 = —1 



2m + 1 



2n+ 1 



(85) 



(86) 
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m 


n 


q/M 


LoM 


oo 


— 


It 


ln3/87r 


-j/8 = 0.04371 


/8 










8 


1 


0.9996 


0.04372 


-i/8 














7 


1 


0.9994 


0.04370 


-i/8 














6 


1 


0.9993 


0.04375 


-i/8 














5 


1 


0.9990 


0.04360 


-i/8 














4 


1 


0.9984 


0.04403 


-i/8 














7 


2 


0.9980 


0.04374 


- 0.12556i 


0.04374 


- 0.37444 










3 


1 


0.9974 


0.04263 


-i/8 














8 


3 


0.9968 


0.04311 


- 0.12364i 


0.04499 


-3V8 


0.04311 


- 0.62636j 






5 


2 


0.9965 


0.04392 


- 0.12334i 


0.04392 


- 0.37666i 










7 


3 


0.9960 


0.04483 


- 0.12348i 


0.04128 


-3i/8 


0.04483 


- 0.62652j 






2 


1 


0.9949 


0.04623 


-i/8 














7 


4 


0.9936 


0.04634 


- 0.12707i 


0.04129 


- 0.37873i 


0.04129 


- 0. 611271 


0.04634 - 


0.87293i 


5 


3 


0.9931 


0.04607 


- 0.12790i 


0.03788 


-3i/8 


0.04607 


- 0. 622901 






3 


2 


0.9919 


0.04484 


- 0.12969i 


0.04484 


- 0.37031i 










4 


3 


0.9903 


0.04222 


- 0.13140i 


0.04860 


-3j/8 


0.04222 


- 0.61860j 






5 


4 


0.9893 


0.04008 


- 0.13200i 


0.04815 


- 0.37830i 


0.04815 


- 0.62170i 


0.04008 - 


0.86800i 


1 


1 


0.9852 


ln2/87r 


- j/8 














3 


4 


0.9780 


0.03849 


- 0.11227i 


0.05121 


- 0.36969i 


0.05121 


- 0.63031i 


0.03849 - 


0.88773i 


2 


3 


0.9743 


0.04193 


- 0.11162i 


0.05331 


-3i/8 


0.04193 


- 0.63838j 






1 


2 


0.9634 


0.04817 


- 0. 111994 


0.04817 


" 0.38801i 










1 


3 


0.9428 


0.05391 


- 0.11419i 


-3i/8 




0.05391 


- 0.63581j 






1 


4 


0.9242 


0.05664 


- 0.11598i 


0.03493 


- 0.35261i 


0.03493 


- 0.64739i 


0.05664 - 


0.88402i 




oo 


Qt 


ln5/87r 


- i/8 = 0.06404 -i 


/8 











TABLE L A sample of roots in cases where the asymptotic spectrum is periodic in the imaginary part. These roots for the 
basis sets out of which the highly damped QNM spectrum is constructed. ^We include the predictions of 16811 for both the 
Schwarzschild limit and the extremely charged black hole limit despite these limits being "singular", cf. the discussion in the 
main text. 

where n and m are integers. One of these roots (for m = 1 and n = 3) is present in Tabled These purely imaginary 
solution may be the cause of some confusion because it is debatable whether they are compatible with the assumptions 
made in the derivation of our QNM condition. In particular, it is not clear to what extent these solutions are relevant 
for a discussion of "purely outgoing wave" solutions. Furthermore, one has to be careful because only some of the 
purely imaginary roots belong to modes that can be traced back to a Schwarzschild QNM. The condition H()8|) was 
derived assuming that the real part of the frequency was positive. Yet, a numerical solution of the condition also 
yields roots with a negative real part. These are not compatible with the underlying assumptions and should be 
discarded as unphysical. These unphysical roots also form spirals as q is varied. Some of the purely imaginary roots 
of the polynomial belong to these presumed unphysical solutions. 

By complementing the polynomial roots with a numerical solution of Ht)8|l . cf. Figure |2| we can further elucidate 
the behaviour of the asymptotic charged black hole QNMs. The numerical solutions illustrated in Figure 13 clearly 
show the spiral nature of the QNMs. Furthermore, we see that the number of times that each spiral touches the 
imaginary axis increases with the "order" of the mode. It is clear that, as Im ljM —>■ — oo the spirals tighten, i.e. the 
QNM frequency becomes exponentially sensitive to variations in q. Whether this is a hint that the problem becomes 
probabilistic in the limit of infinite damping is not clear. 

V. DISCUSSION 

In this paper we have provided WKB results for highly damped quasinormal modes of Schwarzschild and Reissner- 
Nordstrom black holes. In the Schwarzschild case, we have verified previous results for gravitational perturbations. 
In particular, we have provided a clear derivation of the fact that the real part of the quasinormal mode frequencies 
approaches Re ujM = IuS/Stt asymptotically, cf. Eq. (|48|1 . We have also shown that the real part of the asymptotic 
modes vanishes for electromagnetic waves. For Reissner-Nordstrom black holes we have verified the QNM condition 
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recently obtained by Motl and Neitzke f2^, our Eq. H68() . In addition, we have shown that this result holds to leading 
order, not only for the class of perturbations that limits to pure gravitational waves as q —^ 0, but also for the class 
that limits to pure electromagnetic waves. That these two classes of perturbations are isospectral for highly damped 
modes was not known previously. Perhaps the most relevant result concerns the fact that the asymptotic QNM 
frequencies of charged black holes are in general not periodic in the imaginary part. This means that the asymptotic 
spectrum is significantly more complicated than in the Schwarzschild case. We have also explained the breakdown of 
the analysis in the Schwarzschild limit. The WKB approach that we have used should be applicable to many similar 
problems. In particular, it should be possible to use our approach to study the Kerr problem 31, 38]. 

Before we conclude the paper, let us return to the question of a possible link between the QNMs and quantum 
gravity. Do the present results shed further light on this association, or does it now seem as if Hod's original suggestion 
was based on a misleading coincidence? Unfortunately, the available results do not provide a clear answer to this 
question. Hence, we conclude this paper with some speculations. 

First, the original result for gravitational perturbations of Schwarzschild black holes holds. If one accepts the pro- 
posed correspondence between the QNM frequencies and a "transition energy" AM one finds that the quantum of 
area should be AA = 4 In 3. How do results for other perturbing fields fit into this picture? For example, what about 
the fact that the real part of the electromagnetic QNMs vanishes asymptotically? Although interesting in its own 
right, this result is probably not relevant for the discussion. A natural reason for this would be the fundamental dif- 
ference between gravitational and electromagnetic perturbations of Schwarzschild black holes. While electromagnetic 
perturbations correspond to waves propagating in a fixed background, the gravitational waves represent oscillations 
of the spacetime itself. It would perhaps not be surprising if results for the former problem tell us little about the 
quantum levels of the black hole. Could a study of massive fields allow us to make progress? 

The situation is even more complicated when we turn to the Reissner-Nordstrom problem. How can our results be 
understood? At first sight it would be tempting to suggest that the generally nonperiodic nature of the asymptotic 
QNM spectrum, with no unique Re luM , provides an argument against any quantum correspondence. However, this 
conclusion might be premature. One reason for this is that a general perturbation of a Reissner-Nordstrom black 
hole corresponds to a mixture of electromagnetic and gravitational waves. Would it be too surprising if this mixing 
were to prevent a simple correspondence argument? Moreover, we have the unexpected result that we retain the same 
result for the asymptotic QNMs in both the Schwarzschild and the extreme Reissner-Nordstrom limits. Is this telling 
us something profound or is it (again) a mere coincidence? 

Suppose we accept the association between the asymptotic QNMs and the quantum area AA . What would this 
imply for a general charged black hole? The area of a charged black hole is given by 

A = Anrl (87) 

From this we readily get (using n as defined in Eq. (fZ?^ ) 

^ ^ 87rAf(l + k)^ ^ , ^ 87rg(l + k) ^ StttI ^ , ^ 87rgr+ ^ 

A A = ^ —AM — ^Ag = -AM ^—^Aq (88) 

K K ^ Mn Mk ^ ^ ' 

From this relation it is clear that, in general, knowledge of A A/ alone is not sufhcient to "predict" A A. We also need 
to disentangle the electromagnetic waves from the gravitational ones. 

In the absence of more information we can play a simple game. Let us assume that AA is indeed universal and 
thus remains as in the Schwarzschild case (= 4 In 3). Then we can infer from (|88|l that 

. , , kM In 3 , , 

qAq = r+AM 89 

If we associate an "energy" with a given oscillation frequency we have AM — to and thus 

kM In 3 , , 

qAq = r+oj - — 90 

27rr+ 

From this relation we can, given any QNM, compute the value of Aq required for us to obtain the same unique AA. 
In the particular case of extremely charged black holes we would get (since k ^ as g ^ M) 

In 3 

Ag = r+uj = ujM ^ — ^ AM (91) 
87r 

We can also use (|9()|l as a measure of the degree to which a perturbation is "electromagnetic" or "gravitational" . To 
what extent this kind of analysis makes sense is, of course, an open question. 
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Appendix: The Reissner-Nordstrom calculation 

In this Appendix we provide the complete derivation of the WKB condition for highly damped Reissner-Nordstrom 
QNMs. 

The outgoing- wave solution at point a is 

V'a = /^^ (92) 

Moving around the pole at r = r+, cf. Fig. [3 taking full account of all involved Stokes phenomena (assuming that 
the zeros are well isolated) we get 

i^b = /f^' - ift^ = e-^Vf'^ - «e^V^^ (93) 



= e-^'^/f - i [e'^ + e-''^] /^*^^ (94) 

where we have used (|63|l . and then 

= {e*(re+r.) + [1 + e-«(r.+r,)| j-fe) _ ^ ^-2^7] g-^Cr.+r,) ^(t^) 

V'e = {e^(^=+^') + (1 + e^*'') (1 + e-^^^.+r.) | jiu) 

- z (1 e-2n) e-'(^=+^'-^)/^^ (97) 

Moving around the pole at the inner horizon (to e'), the solution does not change because no Stokes lines are 
crossed, but changing the lower limit of integration back to ^3 we get 

V-e' = {e'^= + (1 + e^*^) (1 + e-2n) e-^(r.+2r.) j j(*3) _ ^ ^ ^-2^7) e-'r=/(*3) (gg) 



Then 



= {e^r. ^ ^ g-2»7) [! + (! + e2.7) e-2.r.] g-»r. | ^te) _ ^ _^ ^-2.7) e-»r.^(*3) 
= e-^^ {e^r^ + (1 + e-2n) [l + (l + e^*^) e-^^^^'] 6"*^^} 

- i [l-l-e-2^^] e-*(r<=-^)/i*'^ (99) 



where 



- i(l + e-2n)e-*(^''-^)/i*'' (100) 

V-;, = e-^^ {e*r. ^ ^ g-2.7) [1 + g2^7 + (1 + g2^7) g-2^r,] | ^fe) 

+ ze-*^ {e^r. ^ ^ g-2,7) [! + (! + e^'^) e-'^^'] e'^^^ } /f^ 

= (101) 



A = e^i^^ + (1 + e-2n) (1 + g2.7) (1 + g^2»r,) ^-.r, (^^2) 
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B = 6-2^-^ {e'^= + (1 + [1 + (1 + e^'^) g-^^^,] ^-^T^^ 

= Ae-2''7 _ (1 + e-2n) g-ire (^03) 

Finally returning to the starting point a, we have 

= + (104) 

Thus we see that, in order for the coefficient of /2 to vanish we must require 

A + B = A + _ (1 + = q (105) 

That is, we should have 

A = e-*^= (106) 

which is not too surprising given the result in the Schwarzschild case. As in that problem, one readily shows that the 
(clockwise) monodromy of the solution that is purely ingoing at the event horizon is also e^"^^" . Hence, our calculation 
is consistent. Given the definition of A the WKB condition for highly damped Reissner-Nordstrom QNMs becomes 

e""-^' - 1 - (1 + (1 + e^n) (i + e-2.r.^ (^07) 
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